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Abstract 

The exchange-only optimized effective potential method is implemented with the use of Slater- 
type basis functions, seeking for an alternative to the standard methods of solution with some 
computational advantages. This procedure has been tested in a small group of closed shell atoms 
and diatomic molecules, for which numerical solutions are available. The results obtained with this 
implementation have been compared to the exact numerical solutions and to the results obtained 
when the optimized effective equations are solved using the Gaussian-type basis sets. This Slater- 
type basis approach leads to a more compact expansion space for representing the potential of the 
optimized effective method and to considerable computational savings when compared to both the 
numerical solution and the more traditional one in terms of the Gaussian basis sets. 
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I. INTRODUCTION 



In practical implementations of the Kohn-Sham density functional theory the 
exchange- correlation (XC) energy is usually described by a suitable explicit functional of 
the electron density and the parameters characterizing the density inhomogeneity. A very 
promising step to improve the accuracy of density functional calculations relies on the use 
of orbital-dependent density functionals, in which the XC energy is expressed as an explicit 
functional of Kohn-Sham (KS) orbitals [3]. A rigorous approach to implement the orbital- 
dependent functionals within the KS formalism is the optimized effective potential (OEP) 



method [3H6[, where the XC potential is described by a local multiplicative term and the 
;otal energy functional is orbital-dependent. By virtue of the Hohenberg-Kohn theorems 
l|, the OEP solution is equivalent to the minimization of the total electronic energy with 
respect to the density In the case of the exchange-only — i. e., the Hartree-Fock (HF) 
— energy functional, the corresponding OEP method (xOEP) was first formulated by Sharp 
and Horton Q] and numerically solved in real space for atoms by Talman and Shadwick jsj. 



Applications of real space xOEP formalism to atoms, molecu 
use of pseudopotentials, have been reported in the literature 



es and solids, thanks to the 
8Nl4|. A number of approx- 



imations to the xOEP method, such as the method of Krieger-Li-Iafrate 15|, ll6|, the local 



HF 



17( 1 and the effective local potential 18|] methods have been developed. 



The real space resolution of the xOEP method has only been efficiently applied to hh 



hi 



n 



symmetric systems, such as spherically symmetric atoms and diatomic molecules [5 
Application of the xOEP formalism to polyatomic molecules requires its formulation in terms 
of basis sets suitable for molecular calculations. Currently, there exist several formulations 
of the xOEP method in terms of basis sets of local Gaussian-type (GTO) functions fioj - 



23|. The most popular implementation of the xOEP formalism employs two different basis 



sets, one for the expansion of the KS orbitals and another one for representing the local 
multiplicative potential 22, 2^ 25|. Within this approach, a special care must be taken 
when selecting the auxiliary basis set for the potential, thus leading to the concept of a 
balanced basis set firmly connected to the orbital basis set 25|. Alternatively, a set of the 
products of the occupied and virtual KS orbitals can be employed for the solution of the 



xOEP equations for the local potential 



19 



20 



26|. 



The computational complexity of the xOEP method in a basis set representation depends 
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critically on the size of the orbital expansion basis set. For obtaining faithful solutions of 
the xOEP equations, the orbital basis set should support the linear dependence in the space 
of the occupied- virtual orbital products 



19 



23 



271 . With the use of GTO basis functions, 



this requirement leads to very large orbital basis sets with hundreds of basis functions even 
for small molecules. On the other hand, the use of Slater-type (STO) basis functions give 
considerably more compact orbital basis sets which can be beneficial for the application of 
the xOEP formalism, and a very efficient implementation of quantum chemical formalisms 
with Slater-type basis functions has been achieved in the SMILES suite of programs 



It is the primary purpose of the present work to implement the xOEP method within 
the SMILES package and to analyze the advantages which can be obtained from the use 
of the STO basis sets, testing it for a small number of atoms and diatomic molecules, for 
which both numerical solutions js, 9| and STOs are available 



In this work we employ 



the xOEP algorithm outlined in Refs. 20| and 2l|] and the xOEP equations, formulated 
in terms of the STOs. The solution of the xOEP equation will be carried out through the 
truncated singular value decomposition (TSVD) technique. It will be demonstrated that the 
use of the Slater-type basis sets leads to considerable computational savings in every step 
of the selfconsistent procedure, without deteriorating the accuracy of the calculated xOEP 
total and orbital energies. 



II. THEORY 



In this section the main features of the OEP method will be outlined. In the OEP method 
one seeks for a local multiplicative potential V^(r) such that its eigenfunctions (atomic units 
will be used in this paper) 

|^_iv2 + K(r)^(/)p. = ep,0p, (1) 
minimize the total energy functional given by 

E'^'^'HM] = EE / C(r) (-^V^] M^)dv 

+ J p{r)V,,t{r)dr p(r) J p(r')y^^dTdr' + E,, [{0^4] , (2) 

where runs over occupied orbitals and J2a '^^^^ spin, being p(r) = J2i J2a I0io-P the 
electron density and -Exc [{^i^}] t^e orbital-dependent exchange-correlation energy. Here 
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and below we use indices for occupied orbitals, a,b,... for unoccupied orbitals and 

p,q, ... for general (i.e. occupied or unoccupied) orbitals. In the case of the exchange-only 
formalism, xOEP, the XC energy E^c [{4>ia}] is replaced with the HF exchange energy, 
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The local multiplicative potential 



K(r) = K.t(r) + / T^^dr' + V:{r) (4) 



is splitted into the external potential V"ext(r) (i- e. the potential due to the nuclei), the 
Coulomb potential of the electron cloud and the local exchange potential V^{r). 

The xOEP equations in a basis set representation are obtained from the minimization 



of the total energy presented in Eq. ([2]) with respect to the local potential Vo-(r) [25|, |29| . 
being this minimization equivalent to the minimization respect to the density, by virtue of 
the Sham-Schliiter condition [3] and the Hohenberg-Kohn theorems 1]. 

An expansion of the exchange part of the local potential is assumed in terms of an 



appropiate set of functions [20|, l22|, l2J, l25|, l29|, [30 



V:{r) = J2^,J,{T), (5) 

where a labels the spin, /x labels all needed values of the basis index and w^a are the 
expansion coefficients of the exchange potential in this basis. Following the literature {2^ 
the expansion functions are conveniently defined as 

fU^) = [ f^dr' (6) 



where g^a{'^) are square integrable functions. Note that this definition implies that the 
expansion functions f^aij) are not necessarily square integrable; this is not a problem as the 
local potential does not satisfy this condition. 

The requirement that the total energy to be stationary under the variations of the local 
potential, i. e. 5E^^^^ / 5Vcr{v) = 0, is then equivalent to finding a minimum of the total 
xOEP energy with respect to the set of the expansion coefficients of the the local potential, 
{w^a}- If we work with a real orbital basis and if we introduce a scalar product of two 
functions h and / of our expansion space as 



the minimization of the xOEP energy — i. e. Eg. (|2]) with E^c defined as the Fock exchange 
energy given in Eq. (|3]) — leads to the equation [21 1 



dw 




drdr'4>Ur') [V:{r') - Vr\r,r' 



0. 



Here the non-local potential, V^' (r,r'), is defined as 



(9) 



being {(pqa} the solutions of Eq. ([T]). 
Using the matrix 



in Eq. ([8]) we get a matrix equation equivalent to the minimization of the xOEP energy, 



(10) 



y^xOEP _ 2MM^w - 2Mw"' = 0, 



(11) 



where w is the vector of the expansion coefficients for the local exchange potential and w"' 
is the projection in the chosen basis set of the non-local HP potential, V^'"'''{r, r'). 

In this work our expansion basis set will be a scaled form of the occupied- virtual products, 



specifically 



21 



23 



26| 



9 Mr) 



Thus, the elements of the vector w" are 



w„ 



nl 



drdr 



and the matrix elements M^^^ reduce to 



(12) 



(13) 



ia,jb 



y/^aa ^iay/ ^ba £^ 



(14) 



JO" 



Note that with the basis set of occupied-virtual products it is not possible to get any 
term having a 1/r asymptotic decay. This is corrected by the addition to our basis set 
{g^„{r)} of the Permi-Amaldi function s(r) = p{r)/N, where is the number of electrons 
of the system. That procedure reproduces the Permi-Amaldi potential for long distances and 
will make the xOEP HOMO eigenenergies to be very close, but not equal, to those found 
using the HP method. This prescription {2^ is very different to the one adopted when the 



exchange potential is expanded in an auxiliary basis set 



III. COMPUTATIONAL DETAILS 



As the products of the occupied and virtual states are hnearly dependent, the matrix 
MM^ appearing in Eq. f lTT]) is singular and the equation cannot be solved by inversion 25 |. 
Following the argument given in Refs. 3l| and 27| , only the linearly independent functions 
/^o-(r) (or g^(j(j)) can be used in the expansion of the xOEP, i. e. for a faithful solution of 
the xOEP equation, Eq. f|TT]) . linearly- independent orbital products must be employed. We 
will then apply the Truncated Singular Value decomposition (TSVD) technique to separate 
the linear dependent occupied-virtual products and the independent ones, seeking for the 
linear independent set of products by diagonalization of the MM^ matrix. In order to fulfill 
this condition, a threshold is chosen to discriminate the elements of the MM^ matrix that 
correspond to the linearly independent functions js^]- In this way, the mapping between the 
density and density matrix becomes non-unique and a solution with the energy E^qep > 



27 



32|. 



Ehf is obtained 

In our case, the TSVD method requires the diagonalization of the matrix MM^, in 
general a very large one. But as the STOs represent the unoccupied orbitals in a much more 
efficient way than the GTOs, when Slater-type orbital basis sets are used the size of the 
matrix MM^, and consequently the expansion space for the local potential in Eq. (|5]), is 
much smaller than when using GTOs. This is the main point of this paper: when STO basis 
sets are used, the computational effort (in memory size and in speed of the calculations) for 
each self-consistent cycle is much lower, whereas the quality of the results is preserved. 

The al gori thm outlined in the previous Section was implemented in the SMILES suite of 



programs [28|J, which employs the STO basis sets in quantum chemical calculations. We will 
compare the results of the xOEP calculations obtained with the Slater-type basis functions 
(xOEP-STO) to both numerical exact solutions and xOEP results obtained with the use of 
the Gaussian-type basis functions (xOEP-GTO). The latter results were obtained with the 
use of the MOLPRO2008.1 code 33|], where the xOEP formalism employing Gaussian-type 
basis sets was recently implemented by some of us 20|, |21|. In order to use comparable 
basis sets, the correlation-consistent basis sets of Dunning (cc-pVTZ, cc-pVQZ, and cc- 
pV5Z) were used in the xOEP-GTO calculations, and the STO of similar quality (VBl for 
cc-pVTZ, VB2 for cc-pVQZ and VB3 for cc-pV5Z) [28] were selected for the xOEP-STO 
calculations. These STO and GTO basis sets yield the total HP energies in close agreement 
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(see Tables HI - IIIII below). All the basis sets we have used (xOEP-STO and xOEP-GTO 
calculations) were employed in their uncontracted form. 

Due to the few exact numerical xOEP solutions found in the literature, and to the 
small number of available STOs for atomic and molecular computations, the calculations 
we present here were performed for the Be and Ne atoms and for the LiH, BH, Li2 and CO 
molecules. The numerical solutions were given in Makmal et al. {q] and we have used the 
same internuclear distances (in a.u.): 3.015 for LiH, 2.336 for BH, 5.051 for Li2 and 2.132 
for CO. This distance for the CO molecule has also been used by HeBelmann et al. 22| in 
the xOEP-GTO solution, and we will also compare our results with theirs in Section IIVI 



IV. RESULTS 



The dependence of the xOEP-STO and the xOEP-GTO total energies on the size of the 
basis set and on the TSVD cutoff criterion e for neglecting (near) zero eigenvalues of the 
matrix MM^ has been investigated. The results are collected in Tables HI Hi] and IIIII In the 
calculations, the TSVD cutoff criterion e was varied in the range 10^^ — 10^^. Note that 
when e = 10~^ is employed, the expansion space for the potential is very small because only 
a few eigenvalues of the matrix MM^ are greater than e, yielding energies noticeably above 
the numeric xOEP values. For that reason, we have not reported that energies in the tables. 
Reducing the TSVD cutoff e makes the expansion space bigger and, as a consequence, the 
total xOEP-STO and xOEP-GTO energies decrease, approaching to the accurate numeric 
xOEP values. For e in the range 10~^ — 10~^, the total xOEP energies remain constant to 
within a fraction of mHa. When the xOEP-GTO method is used with a very tight cutoff 
criteria {e < 10~^), the iterative solution of the xOEP equation f lTTj) becomes unstable and 
the xOEP total energy collapses towards the HF energy. As a matter of fact, the procedure 
breaks-down as the matrix MM^ becomes noninvertible. On the other hand, the xOEP- 
STO implementation shows somewhat greater stability and begins to break down at smaller 
values of e, i. e. when e < 10^''. This can be attributed to the fact that, with the use of STO 
functions, the expansion set of the potential is much smaller than with the more traditional 
Gaussian-basis sets (see below). 

The number of the eigenfunctions of the matrix MM^ used for the expansion of the 
potential is also given in Tables [H [TTl and llllt as well as the total dimension of the matrix. 
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It is seen that the dimension of MM^ is considerably smaller in the xOEP-STO method, 
and the dimension of the potential expansion space does not grow as fast as in the case of 
the xOEP-GTO method. So, the xOEP-STO implementation gives a noticeably memory 
savings and a greater stability with respect to the cutoff criterion of the TSVD procedure. 

For all the systems in Tables Ht HTl and UTTt the differences between the xOEP-STO and the 
xOEP-GTO total energies are typically smaller than 1 mHa. When large basis sets are used 
the total xOEP-STO and xOEP-GTO energies approach the exact numeric values with an 
accuracy better than 100 /iHa. It is important to stress here that the xOEP and HP total 
energies converge in a somewhat different way to respect to the basis set size (remember 
that the basis sets we use in this paper were nor developed neither optimized for the xOEP 
calculations) so the differences E^qep — Ehf oscillate. 

Table [IV] summarizes the results of calculations for Be and Ne atoms and a number of 
diatomic molecules studied by Makmal et al. [9] using the real space xOEP method. For 
the sake of comparison, the xOEP energies obtained by HeBelmann et al. {22] with the use 
of balanced auxiliary basis sets for the potential are also shown when available. Note that 
the xOEP-STO energies are typically in a somewhat better agreement with the numerical 
values than the xOEP-GTO energies obtained with similar basis set (the STO basis sets 
give results about 0.1 mHa below the energies obtained by HeBelmann et al.). 

Table |V] collects the energies of the occupied orbitals obtained with the use of the xOEP- 
STO and xOEP-GTO methods for the Be atoms and the LiH and Li2 molecules. There is 
a good agreement between the numerical values of these orbital energies and the energies 
obtained with the two xOEP methods. Furthermore, the orbital energies from the xOEP- 
STO and from the xOEP-GTO calculations agree with each other to within a few mHa. 

Besides the total and orbital xOEP energies, we have also studied the energy decom- 
position into the kinetic, nuclear-electron attraction, electron-electron repulsion and the 
exchange energies. Table I VI I presents the above components of the total xOEP energy as 
obtained using the STO and GTO basis sets for the Be atom and the CO molecule (there 
are no numerical solutions available in this case for the xOEP). 

Figs. [1] and |2] show the results of the xOEP potential for the Ne atom and the CO molecule 
along the main axis. For the Ne atom, note the close agreement between our result and that 
obtained with the exact numerical calculation by Kurth and Pittalis 3J|, thus yielding a 
smooth potential that shows only small deviations from the numerical potential; in any 
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case, these xOEP potentials are also very similar to that evaluated with a Gaussian basis 



set [31j using the procedure presented in Refs. [20| and |2l|. In the lower panel of Fig. [2] the 
results we have obtained using the VB3 basis set for the CO molecule are compared with 
the calculation by HeBelmann et al. [22|, using an auxiliary basis set within a Gaussian 
representation; the STO results show a good agreement with the xOEP-GTO potential. 
For the sake of completeness, in the upper panel we also present several other calculations 
for the internuclear region, using STO basis sets of different quality (VBl, VB2 and VB3). 
Note that our xOEP potentials do not present any unphysical wiggle as those found by 



Staronerov et al. 



18] and have a good agreement with both the numerical calculation and 



the xOEP-GTO solution by Hessellmann et al. [22| with the use of auxiliary basis set. 

In summary, the previous results show that both xOEP-STO and xOEP-GTO methods 
yield close results to the numerical exact ones for the total energies, the one-electron energies 
of the occupied orbitals and the xOEP potentials. 

V. CONCLUSIONS AND DISCUSSION 

The implementation of the xOEP formalism with Slater-type basis functions has been de- 
veloped. This procedure has been tested in a small group of closed shell atoms and diatomic 
molecules, for which both numerical xOEP solutions and STO basis sets are available. When 
compared to the exact numerical solution of the xOEP equations we have obtained very good 
results; they are even a bit better than those given by the Gaussian-type basis sets procedure. 

O and xOEP-GTO results obtained with the prescrip- 



On the other hand, both xOEP-S 

n r 

tion proposed in Refs. [20| and 2l| give energies about 0.1 mHa below those obtained for 



22| 



xOEP-GTO by HeBelmann et al., and thus they are closer to the exact results 

The new method leads to a considerably more compact expansion space for the xOEP 
local multiplicative potential, yielding noticeable savings in the computational effort to be 
done in each one of the cycles of the self-consistent procedure. Yet another advantage of 
using the Slate- type basis sets is that, within the TSVD algorithm, fewer eigenvalues of 
the (near) singular matrix MM^ need to be employed, which leads to an increased numeric 
stability of the xOEP-STO method as compared to the xOEP-GTO algorithm. 

As a final remark, it is known that a more efficient xOEP algorithm can be developed 
based on the use of the incomplete Cholesky decomposition technique. The use of this 



9 



technique would facilitate the apphcation of the xOEP-STO method to larger molecules. 
This implementation is currently in progress and will be reported elsewhere. 

On the other hand, this work is a first step to develop a local potential formalism for 

both exchange and correlation. Due to the smaller memory requirements of the STO scheme 
we have presented here, it can be used to study bigger molecules than those that can be 
solved with the standard GTO approach. The implementation of the correlation part of the 
potential is under development. 
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TABLE I: Total energies (Hartree) for the Be atom using the HF and the xOEP methods with the 
STO and the GTO basis sets. In the first column we indicate the value e used as a threshold in the 
TSVD decomposition of the TSVD decomposition of the MM^ matrix in the xOEP method. For 
the STO and GTO xOEP results we have included, within brackets, the total number of occupied- 
virtual products (first number in the brackets) and the number of them that are used (second 
number in the brackets) in each calculation. 



e 


STO~basis 


jpSTO 


GTO-basis 






VBl 


-14.57238 (12/3) 


cc-pVTZ 


-14.57233 (26/10) 




VB2 


-14.57240 (17/4) 


cc-pVQZ 


-14.57234 (38/17) 




VB3 


-14.57240 (22/5) 


cc-pV5Z 


-14.57227 (52/20) 


10-^ 


VBl 


-14.57242 (12/4) 


cc-pVTZ 


-14.57233 (26/14) 




VB2 


-u.f)nm (17/5) 


cc-pVQZ 


-14.57256 (38/20) 




VB3 


-14.57245 (22/7) 


cc-pV5Z 


-14.57255 (52/26) 




VBl 


-14.57244 (12/6) 


cc-pVTZ 


-14.57241 (26/16) 




VB2 


-14.57256 (17/8) 


cc-pVQZ 


-14.57256 (38/21) 




VB3 


-14.57254 (22/13) 


cc-pV5Z 


-14.57255 (52/30) 




STO-basis 




GTO-basis 






VBl 


-14.57297 


cc-pVTZ 


-14.57287 


HF 


VB2 


-14.57298 


cc-pVQZ 


-14.57296 




VB3 


-14.57301 


cc-pV5Z 


-14.57301 



variational energy obtainable with the functional given in Eq. ([3]), that is the HF energy Ehf- 
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TABLE II: Total energies (Hartree) for the LiH molecule using the HF and the xOEP methods with 
the STO and the GTO basis sets. In the first column we indicate the value e used as a threshold 
in the TSVD decomposition of the MM^ matrix in the xOEP method. For the xOEP-STO 
and xOEP-GTO results we have included, within brackets, the total number of occupied-virtual 
products (first number within the brackets) and the number of them that are used (second number 
within the brackets) in each calculation. 



e 


STO-basis 


rpSTO 


GTO-basis 


rpGTO 


10-3 


VBl 


-7.98651 (14/3) 


cc-pVTZ 


-7.98660 (58/24) 




VB2 


-7.98660 (18/4) 


cc-pVQZ 


-7.98661 (88/32) 




VB3 


-7.98669 (26/6) 


cc-pV5Z 


-7.98671 (134/42) 


lo-'i 


VBl 


-7.98658 (14/6) 


cc-pVTZ 


-7.98661 (58/29) 




VB2 


-7.98693 (18/7) 


cc-pVQZ 


-7.98698 (88/44) 




VB3 


-7.98G93 (26/9) 


cc-l)V5Z 


-7.98696 (134/63) 


10-5 


VBl 


-7.98681 (14/8) 


cc-pVTZ 


-7.98678 (58/36) 




VB2 


-7.98699 (18/9) 


cc-pVQZ 


-7.98701 (88/52) 




VB3 


-7.98689 (26/13) 


cc-pV5Z 


-7.98700 (134/74) 




STO l)asis 


T7STO 


GTO l)Msis 


jtGTO 




VBl 


-7.98479 


cc-pVTZ 


-7.98695 


HF 


VB2 


-7.98650 


cc-pVQZ 


-7.98722 




VB3 


-7.98726 


cc-pV5Z 


-7.98733 
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TABLE III: Total energies (Hartree) for the CO molecule using the HF and the xOEP methods 
using the STO and the GTO basis sets. In the first column we indicate the value e used as 
a threshold in the TSVD decomposition of the TSVD decomposition of the MM^ matrix in the 
xOEP method. For the xOEP-STO and the xOEP-GTO results we have included, within brackets, 
the total number of occupied-virtual products (first number within the brackets) and the number 
of them that are used (second number within the brackets) in each calculation. 



£ 


STO-basis 


rpSTO 


GTO-basis 


TpGTO 


10-^ 


VBl 


-112.77681 (93/10) 


cc-pVTZ 


-112.77573 (438/94) 




VB2 


-112.78251 (111/13) 


cc-pVQZ 


-112.78246 (686/104) 




VB3 


-112.78433 (184/21) 


cc-pV5Z 


-112.78428 (1046/116) 


10-4 


VBl 


-112.77999 (93/16) 


cc-pVTZ 


-112.77657 (438/128) 




VB2 


-112.78469 (111/20) 


cc-pVQZ 


-112.78392 (686/176) 




vb;3 


-112.78521 (181/25) 


cc-pV5Z 


-112.78519 (1016/238) 


10-5 


VBl 


-112.77780 (93/24) 


cc-pVTZ 


-112.77704 (438/160) 




VB2 


-112.78524 (111/27) 


cc-pVQZ 


-112.78464 (686/226) 




VB3 


-112.78526 (184/32) 


cc-pV5Z 


-112.78551 (1046/312) 




STO l)asis 


T7STO 


GTO l)asis 


rGTO 




VBl 


-112.78199 


cc-pVTZ 


-112.78014 


HF 


VB2 


-112.78959 


cc-pVQZ 


-112.78891 




VB3 


-112.79056 


cc-pV5Z 


-112.79064 
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TABLE IV: Total xOEP energies (Hartree) for several atomic and diatomic systems calculated 
using a TSVD decomposition for the TSVD decomposition of the MM^ matrix in the xOEP 
method (the threshold e used to discriminate the eigenvalues was 10~^). 





Numeric 


ST0(VB3) 


GT0(cc-pV5Z) GTO(Hess)'' 


Be 


-14.57254'' 


-14.57254 


-14.57255 


-14.57243 


Ne 


-128.54553'^ 


-128.54540 


-128.54548 


-128.54538 


LiH 


-7.98691^ 


-7.98689 


-7.98700 




Li2 


-14.87076" 


-14.87044 


-14.87090 




BH 


-25.12963'^ 


-25.12941 


-25.13013 




CO 


-112.785(3)" 


-112.78526 


-112.78551 


-112.78491 



"Ref. 
^Ref. 
'^Ref. 
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TABLE V: Orbital energies (a. u.) for the Be atom and the LiH and Li2 molecules. Numerical 
results for the xOEP are presented. All STO (GTO) calculations have been done with the VB3 
(cc-pV5Z) basis set. The threshold e used in the TSVD to discriminate the eigenvalues of the 
MM"!" matrix in the xOEP method was 10 . 





HF-STO HF-GTO 


xOEP-numeric 


xOEP-STO xOEP-GTO 


Be 






Is 


-4.17045 


-4.17326 


-4.1668 


-4.1689 


-4.1704 


2s 


-0.34903 


-0.34925 


-0.34885 


-0.34878 


-0.34892 


LiH 






la 


-2.44543 


-2.44534 


-2.0786 


-2.09092 


-2.07071 


2a 


-0.30172 


-0.30172 


-0.3011 


-0.31391 


-0.31384 


Li2 






lag 


-2.45311 


-2.44994 


-2.0276 


-2.00774 


-2.01361 


ICTu 


-2.45311 


-2.44994 


-2.0272 


-2.01262 


-2.00738 


2ag 


-0.18194 


-0.18193 


-0.1813 


-0.18401 


-0.18616 
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TABLE VI: Decomposition of the total energy (Hartree) in its terms for the Be atom and the CO 
molecule. All STO (GTO) calculations have been done with the VB3 (cc-pV5Z) basis set. The 
threshold e used in the TSVD to discriminate the eigenvalues of the MM^ matrix in the xOEP 
method was 10~^. 





HF-Numeric 


HF-STO xOEP-STO 


HF-GTO xOEP-GTO 


Be 








Ts 
E 

TP a 
Etot 


-14.57304^ 


14.57301 14.57319 
-33.63509 -33.63271 

4.48911 4.48964 
-14.57300 -14.57254 


14.57301 14.57309 
-33.63518 -33.63387 

4.48916 4.48944 
-14.57301 -14.57255 


CO 








Ts 
E 

TP a 

Etot 


-112.79078'= 


112.64045 112.63314 
-310.87673 -310.86921 
62.93165 62.93421 
-112.79056 -112.78526 


112.64200 112.63348 
-310.88020 -310.86944 
62.93349 62.93637 
-112.79064 -112.78551 



'^This 
''Ref. 
'^Ref. 



teltm is the sum of the Hartree and the exchange energies. 
36| 
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0.01 0.1 1 10 

r(a.u.) 

FIG. 1: xOEP potential for the Ne atom. Atomic units are used in both coordinate axis as well as a 
logarithmic scale in the OX axis. The dashed blue line corresponds to the xOEP-GTO calculation 



of Fernandez et al. 



31[, the red line represents the results of this paper and the green one reflects 



the numerical results of Kurth and Pittalis 



34l |. In the inset (where the OX scale is now linear) 



we show in more detail the region around the shoulder. 
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FIG. 2: xOEP potential for the CO molecule along the axis of the molecule. Atomic units are 
used in both coordinate axis (we have used an internuclear distance of 2.132 a.u.). In the lower 
panel, the red line show the results of this paper and the green one reflects the xOEP-GTO results 
of HeBelmann et al. \^\. In the upper panel the region between the nuclei of the molecule is 
depicted in more detail; the xOEP-GTO results are shown again in green, and results for xOEP- 
STO calculations with different STO basis sets [28|] are also plotted (VBl, blue dashed line; VB2, 




magenta dotted line; VB3, red thick line). 
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